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Abstract 

Quark deconfinement phase transition at finite temperature and 
density is investigated in the frame of quantum mechanics. By solv- 
ing the Schrodinger equation for a heavy quark in a thermal mean 
field, we calculate the quark probability distribution as a function of 
temperature and density. The confined wave function in the vacuum 
expands outward rapidly when the temperature and density are high 
enough. The obtained phase transition line agrees qualitatively with 
the result of lattice QCD. 

It is widely accepted that the collective effect of a multiparton system can 
change the vacuum structure of Quantum Chromodynamics (QCD)pQ, and 
the partons confined in a hadron bag can move out when the temperature 
and density are high enough and form a new state of matter, the so called 
quark-gluon plasma (QGP)[2j. This kind of new matter may be produced in 
the early stage of a high-energy heavy-ion collisionjH]. The dense partons, 
produced in the vacuum at high temperature or compressed at high baryon 
density, make the quark wave functions overlap in a space with dimension of 
the colliding nuclei. An important question in relativistic heavy ion collisions 
is how to identify the new state of matter if it is created in the early stage. 
According to the numerical calculation of lattice QCD [3], the critical tem- 
perature of quark deconfinement is about 150 MeV when the density effect 
is excluded. 

The idea that the overlap of wave functions, or the space extension of 
wave functions at finite temperature and density is the original reason of 
quark deconfinement phase transition arises naturally an interesting ques- 



tion to study the behavior of quark distribution function in the frame of 
quantum mechanics. In this letter, we will study a quark moving in a ther- 
mal mean field describing multiparton interactions, and obtain the quark 
probability distribution by solving the Schrodinger equation at finite tem- 
perature and density. The temperature and density at which the quark wave 
function expands outward rapidly can be considered as the critical point of 
quark deconfinement. Because of the limitation of non-relativistic quantum 
mechanics, we focus on the motion of a heavy quark at finite temperature 
and density, corresponding to, for example, the dissociation of J/ip[5\ as a 
signature of QGP in relativistic heavy ion collisions. 

In the string-like models[3 that describe the quark confinement well, 
one uses a linear potential V — kr to express the interaction between two 
quarks. To simplify the numerical calculation in the following, we use a three 
dimensional square well 



instead of the linear potential. Here a is the boundary of the region where 
the motion of the quark is asymptotically free, and can be taken as the 
radius of the hadron constructed by the considered quarks, b is the space 
scale of the system, and Vq should be so large that the quark wave function 
can be confined in the asymptotically free region r, a. Since the region of 
collective interaction in high-energy heavy-ion collisions is of the order of the 
colliding nuclei, we treat b as the radius of a heavy nucleus. In our numerical 
calculations, we take the parameters of the confinement potential as a = 1 
fm, b = 5 fm, and V = 10 GeV. 

At finite temperature and density, the quark will be affected by many- 
body interactions in the region r > a in addition to the confinement po- 
tential. In mean field approximation, the thermodynamic potential of the 
quasi-particle system is [I] 
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with the quark and antiquark densities 

f q {T,n) = , 

e t +1 

MT,ii) = , (3) 
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where g is the quark degree of freedoms of flavors, colors and spins, T and 
ji are respectively the temperature and baryon chemical potential, E q = 
yj miff + V 2 is the quark energy with effective mass m e ff obtained in the 
mean field. Because the contribution from heavy quarks to the thermody- 
namic potential is very small compared with light quarks, we will only take 
u and d quarks into account, i.e., g = 12. 

From the rules of thermodynamics, the pressure of the system is 

P(7» = -fi(7» . (4) 

— * 

Considering the relation between the potential V th and the force F th that the 
heavy quark we study feels in the region r > a, 

F a (r) = -Wo,(f) , (5) 

-^ = -4^P(T,„), (6) 

and assuming that the potential falls down with increasing r and vanishes 
on the boundary of the system, namely V t h(r = b) = 0, the potential V t h can 
be expressed as 

V th (r\T, v) = -f r b dr'4nr' 2 P(T, //) = ^tt(6 3 - r 3 )Q(T, /i) . (7) 

With the obtained potential, the motion of the heavy quark at finite tem- 
perature and density is described by the Schrodinger equation 

i-^mt) = [-—V 2 + V ef f(r\T, f i))^f,t) , (8) 

where m is the mass of the heavy quark, and V e ff is the effective potential 
defined as 

[ < r < a 

V e f f (r\T,fi) = I V + V th (r\T,fi) a<r<b (9) 
[ oo r > b 
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By separating the variables of the wave function ip(r, t) = R(r)Y(6, <p)e tEt , 
and letting <f>(r) = rR(r), we obtain the radial equation for a stationary state 

1 ^ W + ( V e „(r\T, „) + #M = E4,(r) , (10) 
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where E is the eigenvalue of the quark energy. Since we focus in this letter on 
quark deconfinement at finite temperature and density, we discuss in the fol- 
lowing the ground state of the Schrodinger equation only. The corresponding 
radial equation is then simplified as 

1 d 2 Mr) 



2m dr 2 



+ V eff (r\T,(jL)Mr)=E Mr) • (H) 



We first discuss qualitatively the deconfinement phase transition through 
the behavior of the effective potential as a function of T and fi. From 
the T- and /i-dependence of the thermodynamic potential, we can imme- 
diately obtain the conclusion that the confinement potential in the vacuum 
at T = n = will be suppressed by the thermal mean field V t h < 0, and the 
motion of the quark will extend outside when the temperature and density 
increase. Since the mean field decreases monotonously from the maximum 
on the boundary r = a to the minimum on the other boundary r = b, the 
suppression at r = a is most strong, and the potential at r = b does not 
change at finite temperature and density. 

In order to illustrate quantitatively the above estimated quark deconfine- 
ment phase transition, we should solve the Schrodinger equation (jll)) at finite 
temperature and density to get the quark probability distribution 

p(r\T,fx) ocr 2 R 2 (r\T,fx) . (12) 

We take the considered heavy quark in the Schrodinger equation as a c quark 
which forms J ftp with m = 1.6 GeV, and the effective mass of the light u 
and d quarks in V e ff as their constituent quark mass m e ff = 0.35 GeV. 

In the limit of fi = 0, the probability distribution as a function of tem- 
perature obtained by solving numerically the Schrodinger equation (|TT|) is 
shown in Fig.(£Q). In the vacuum with T = \x = 0, the chosen parameters for 
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the square well guarantee that almost the whole quark wave function is re- 
stricted in the region r < a, that means quark confinement. As temperature 
increases the quark motion extends outside very slowly in the beginning, then 
the extension is accelerated rapidly in the neighborhood of T = 120 MeV, 
which indicates that the two c quarks forming J / ip begin to dissociate in the 
hot mean field, and finally the quark moves in the whole region of the system 
< r < b. The case T = is similar to that of ii = 0, the quark probability 
distribution as a function of r for different chemical potentials is indicated 
in Fig. (|2*]) . The probability expands outward rapidly in the neighborhood of 
ll = 445 MeV. 
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Figure 1: The quark probability distribution at finite temperature for fj, = 0. 
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Figure 2: The quark probability distribution at finite chemical potential for 
T = 0. 



It is well known that for a particle moving in a central field the most 



probable radius describes the bound degree of the stationary state. The 
most probable radius r (T, /i) of the c quark motion at finite temperature 
and density defined by 

dp(r\T, n) 



dr 



\r=r 







(13) 



is indicated as a function of T for fi = in Fig.(j3J). It is clear to see that 
around T = 120 MeV tq goes up from a/2 to 6/2 very fast. We can de- 
termine the critical temperature T c and chemical potential fi c of the quark 
deconfinement phase transition as the point at which the most probable ra- 
dius r (T c , /i c ) reaches a, the boundary of the asymptotically free region in 
the vacuum, 

r (T c ,/i c ) = a . (14) 

The critical line which separates the deconfinement phase from the confine- 
ment phase is shown in Fig.(0J). While the most probable radius ro jumps 
from a/2 to 6/2 in a very narrow region around T c , it is still continuous at the 
critical point, and the phase transition we study is a continuous one. When 
the temperature and chemical potential exceed the critical values, the region 
where the quark is asymptotically free will expand, i.e., the position of a will 
move outward. 
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Figure 3: The most probable radius as a function of temperature for /i = 0. 



In summary, by describing the collective effect as a thermal mean field 
and combining it with the confinement potential in the vacuum, we have 
demonstrated the physical picture of quark deconfinement in the frame of 
quantum mechanics through solving the Schrodinger equation and calculating 
the quark probability distribution at finite temperature and density. When 
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Figure 4: The phase transition line separating the region with quark confine- 
ment from that with quark deconfinement. 

the temperature and density are high enough, the quark wave function ex- 
tends outward from inside the hadron. The critical temperature and density 
around which the wave function extension happens rapidly separates the con- 
finement from deconfinement. The obtained critical temperature T c = 120 
MeV at n — agrees qualitatively with the lattice QCD calculation. 

Acknowledgments: The work was supported by Grant Nos. NSFC19925519, 
10135030, and G2000077407. 

References 

[1] J. I. Kapusta, Finite Temperature Field Theory, Cambridge University 
Press, 1993; M. le Bellac, Thermal Field Theory, Cambridge University 
Press, 2000; 

[2] Quark-Glon Plasma, ed. R.C.Hwa (World Scientific, 1990). 

[3] C. Y. Wong, Introduction to High-Energy Heavy-Ion Collisions, World 
Scicetific, 1995. 

[4] F. Karsch, Nucl. Phys. A698(2002)199c. 

[5] T. Matsui and H. Satz, Phys. Lett. B178 (1986) 416; C. Gerschel and 
J. Hiifner, Ann. Rev. Nucl. Part. Sci. 49 (1999) 255. 



7 



